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Abstract 

We present a study, within a mean-field approach, of the kinetics of a 
classical mixed Ising ferrimagnetic model on a square lattice, in which the two 
interpenetrating square sublattices have spins a = ±1/2 and S = ±1,0. The 
kinetics is described by a Glauber-type stochastic dynamics in the presence 
of a time-dependent oscillating external field and a crystal field interaction. 
We can identify two types of solutions: a symmetric one, where the total 
magnetization, M, oscillates around zero, and an antisymmetric one where 
M oscillates around a finite value different from zero. There are regions of 
the phase space where both solutions coexist. The dynamical transition from 
one regime to the other can be of first or second order depending on the 
region in the phase diagram. Depending on the value of the crystal field we 
found up to two dynamical tricritical points where the transition changes from 
continuous to discontinuous. Also, we perform a similar study on the Blume- 
Capel (5 = ±1,0) model and found strong differences between its behavior 
and the one of the mixed model. 



I. Introduction 

The time evolution of metaestable states of a mixed Ising ferrimagnetic model are 
studied by establishing a Glauber-type dynamic ||l[ that drives the system between 
two equivalent ordered phases. This approach provides a simple way to introduce 
a nonequilibrium kinetics into models such as the Ising and the Ising mixed sys- 
tems, that do not have a deterministic dynamics. Numerous studies indicate that 
stochastic dynamics, despite it simplicity, can describe in a qualitative correct way 
nonequilibrium phenomena and dynamical phase transitions found in real systems 
1^. Recently a kinetic Ising model has been applied sucesfuUy to model magnetiza- 
tion switching in nanoscale ferromagnets 
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The model to be studied is a mixed Ising ferrimagnetic system on a square 
lattice in which the two interpenetrating sublattices have spins one-half (±1/2) and 
spins one (±1,0). This system is relevant for understanding bimetallic molecular 
ferrimagnets that are currently being synthesized by several experimental groups in 
search of stable, crystalline materials, with spontaneus magnetic moments at room 
temperature 0. The magnetic properties of the mixed Ising models have been 
studied by high-temperatures series expansions P], renormalization-group 0, mean- 
field 1^, Monte Carlo simulations and numerical transfer matrix calculations 0, but 
always at equilibrium conditions. 

In this work we are going to study, within a mean field approach, the kinetics of 
de model in presence of a time-dependent oscillating external field. We found that 
there are three regimes: a paramagnetic one were the total magnetization follows 
the field, a ferromagnetic one where the total magnetization does not follows the 
field and oscillates around a value different from zero, and a region where both so- 
lutions coexist. The transition between regimes can be continuous or discontinuous, 
depending on the temperature and the external field. Depending on the parameters 
of the Hamiltonian this system can have up to two dynamical tricritical points. 

II. The Model 

We consider a kinetic mixed Ising ferrimagnetic system in a square lattice de- 
scribed by the Hamiltonian 

H = -J E S^cr, -Dj2s^- HiJ2S^ ± E^.) (1) 

(nn) i i j 

where the Si take the values ±1 or and are located in alternating sites with spins 
(Tj = ±1/2. Each spin S has only cr-spins as nearest-neighbors and vice- versa. The 
sum I]<nn> is carried out over all nearest-neighbors pairs. The sums X^i and 
run over all the sites of the S and a sublattices, respectively. J is the exchange 
parameter, D is the crystal field interaction, and if is a time-dependent external 
magnetic field given by 

H(t) = Hocos{ujt) , (2) 

all in energy units. We choose J to be negative, so the coupling between the nn 
spins is antiferromagnetic. 

The system evolves according to a Glauber-type stochastic dynamics |]1|, that 
is described a continuation, at a rate of 1/r transitions per unit time. Leaving 
the S'-spins fixed, we define P'{ai, . . . ,aN',t) as the probability that the system 
has the cr-spins configuration, ai, . . . , ctat, at time t, also, by leaving the cx-spins 
fixed, we define P"{Si, . . . , Sjy; t) as the probability that the system has the S'-spins 
configuration. Si, ... , Sn, at time t. Then, we calculate Wj{aj) and Wi{Si S[), 
the probabilities per unit time that the jth a-spin changes from aj to —aj and the 
ith S'-spin changes from Si to S'^, respectively. 

II. a. Calculation of Wj{aj) 
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We write the time derivative of P'{<Ji, ■ ■ ■ , ctn] t) as 

^P'(ai, . . . , (Tat; t) = -[^ Wj{crj)]P'{ai, . . . , aj, . . . , aN] t) 

+ J2 Wji-(^j)P'i(^u ■■■,-crj,...,crN;t) . (3) 

j 

If the system is in contact with a heat bath at temperature T, each spin a can 
flip with probabihty per unit time given by the Boltzmann factor 

1 exp(-/3AEj) 
"-•^"^•^ = r l + exp(-/3Ai^-) ' /^=Vfc-^ 

where 

A£;; = 2a,(J^5, + /7) (5) 

<i> 

gives the change in the energy of the system when the aj spin flips. 

From the master equation associated to the stochastic process it is simple to 
prove that the average < crj{t) > satisfies the following equation |10| 

< (Tj >= - < a, > + < i tanh[^/?( J Si + H)] > , (6) 

that, within a mean-field approach |ll| and for an external field defined by (|^), takes 
the form 

d 11 

T— < a >= - < a > +- tanh ^(3[JZ < S > +Ho cos(u;t)] , (7) 

where < a > and < S > are the a-spin and S-spin sublattice magnetizations 
respectively, and Z = 4 is the coordination number for this model. 

Il.b. Calculation of Wi{Si Si) 

In a similar way, the time derivative of the P"{Si, . . . , SN',t) can be written as 

^P"(5i, . . . , t) = - E[ E W.{S. S:)]P"{S,, . . . , 5^; t) 

+ E[E W.{S:-^S,)P"{Si,...,S:,...,SN;t)]. (8) 

Each spin can change from the value Si to the value SI with probability per unit 
time 

Wfc, c.^ _ 1 exp{-f3AE"{S. ^ Sj)) 

vvi[bi - -— — — [y 

i"Es^'exp(-/5AE"(S^i S'i)) 

where the J2s' is the sum over the 3 possible values of S'^, ±1,0, and 

AE"{S, ^ Si) = -{S[ - S.){J E a, +H)- (^'J - S^)D (10) 

<j> 
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gives the change in the energy of the system when the Si spin changes. The proba- 
bihties satisfy the detailed balance condition 



Wi{Si SI) _ Peq{Sl, . . . , S*-, . . . , Sn) 



^ii^'i ~^ ^i) Peq{^l^ ■ ■ ■ , Si, . . . , S^ 



N) 



and, substituting the possible values of Si, we get, 

exp{-f3D) 



W,{1- 


- 0) 


= W,{-1- 


-0) 


W,{1- 


^-1) 


= WiiO ^ 


-1) 


WiiO- 


- 1) 


= W,{-1- 


-1) 



2 cosh a + exp( 




exp(— a) 




2 cosh a + exp( 


-m 


exp(a) 




2 cosh a + exp( 


-m 



(12) 



where a = P{JJ2<j> c"j + H). Notice that, since Wi{Si — >• S^) does not depend on 
Si we can write, Wi{Si S"-) = Wi{S[), and the master equation becomes 

^P"iS,, . . . , 5^; t) = - E[ E W.iS'^]P"iSu ...,S^;t) 

+ T.W^iS^)[Y. P"{Su...,S[,...,SN;t)] (13) 

Since the sum of probabilities is normalized to one, by multiplying both sides of 
(|l^) by Sk and taking the average, we obtain 

r A < 9 >= - < 9 > I 2sinh/3(JE<,>^. + g) . . 

dt ' ' 2cosh/3(JE<,>cT, +i/) +exp(-/3Z}) ^ ' 

or, in terms of a mean-field approach 

d ^ ^ 2smh/3\JZ < a > +Hocos(iut)] 

dt 2 cosh /3[JZ <a> +Hq cos(cjt)] + exp(-/5L)) ' ^ ' 

III. Results 

The system evolves according to the set of coupled differential equations given 
by (0) and (|l|), that can be written in the following form. 



^—mi = -mi + sinh^(m2 + /icosQ 

d^ coshj;{m2 + hcos^) + ^exp{—/3D) 

d IX 
^-^^2 = -1712 + - tanh — (mi + h cos ^) 

where mi =< 5* > and 1712 =< cr >, Q = tuj, ^ = ujt, T = {[3JZ)~^ and h = Hq/JZ. 
We fix J = —1. We are going to study the stationary solutions of this system and 
its dependence with the parameters. 
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In general, the solutions of the system depend on the initial conditions, = 0) 
and m2(^ = 0). For a given set of parameters and initial conditions the system pass 
through a transient regime until the solution becomes stationary. In Fig. 1, we 
show the solutions of the total magnetization, M = mi + m2, for different initial 
conditions. We found that the following combinations of initial conditions, 

mi(^ = 0) = l, m2(e = 0) = 0.5 (17) 
mi(e = 0) = 0, m2(e = 0)=0 (18) 

give all the possible different stationary solutions. As expected, the sublattice mag- 
netizations, mi, 7712 and the total magnetization M are periodic functions of ^ with 
period 27r. 

By choosing as the dinamical parameter the average total magnetization in a 
period, 

1 



Q = — MiOd^ (19) 

ZTT JO 

we can identify two types of solutions: a symmetric one where M(^) follows the field, 
oscillating around zero giving Q = 0, and an antisymmetric one where M(^) does not 
follows the field and oscillates around a finite value different from zero, such that Q ^ 
0. Examples of both types of solutions are shown in Fig. 2. The first type of solution 
is called paramagnetic (P) and the second type ferromagnetic (F). The dynamical 
phase transition is located in the boundary between both solutions, and depending 
on the region of the phase diagram, can be of first or second order. For the first order 
phase transition the order parameter Q is discontinuous, and jumps abruptly from 
zero to a nonzero value. For the second order phase transition, the order parameter 
decreases continuously to zero, but its first derivative is discontinuous. If we fix the 
values of D and 1] for each value of we get that the dynamical phase transition 
occurs at temperature T^. In the thermodynamical plane, {T,h), we can define a 
critical line = Tc{h), the point on the critical line at which the transition changes 
from first order to second order is denominated dynamical tricritical point. 

As it is shown in Fig. 3, for D large and positive, we recover the phase diagram of 



the standard Ising model |Tl|] . At high temperatures the solutions are paramagnetic 
and at low temperatures there are ferromagnetic. The boundary between both 
regions, (F)— i>(P), is given by the critical line, and indicates a continuous phase 
transition. At low temperatures, there is a range of values of h for which the 
ferromagnetic and paramagnetic solutions coexist and two separated critical lines 
appear, one that indicates the discontinuous transition between the (F) and the 
(P+F) regions, and the other that indicates the discontinuous transition between 
the (P+F) and the (P) regions. The point where both lines merge signals the change 
from a first order to a second order phase transition and is a dynamical tricritical 
point. 

When D becomes negative, we found that in the range —2.03760 < D < 
—0.11825 (for 1^ = 1) and low temperatures, a new first order critical line ap- 
pears that separates a new region where both solutions coexist (F+P) from the (F) 
region. A phase diagram corresponding to a value of D in this interval is shown 
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in Fig. 4. This new critical line is due to the fact that, in this region, the sub- 
lattice magnetization mi is very small, almost zero, such that both sublattices are 
practically uncoupled, leaving m2 free to oscillate with the external field giving a 
paramagnetic type response, Q = 0, as seen in Fig. 5. When D < —2.03760 the 
phase diagram changes dramatically. The paramagnetic region is extended over the 
low h regions, new critical lines appear and a second dynamical tricritical point 
emerges. Phase diagrams in this interval can be observed in Fig. 6. The two first 
order transitions that occur at low temperature and the two second order transitions 
at higher temperature are shown in Fig. 7. In all cases we found that, as the crys- 
talline field D becomes larger but negative, the region of the phase space where the 
system behaves ferromagnetically becomes narrower and shifts toward higher values 
of h. This effect, shown in the inset of Fig. 6 is somewhat expected since, for these 
values of D, higher values of h are needed to keep mi(^) different from zero such 
that the order parameter, Q, can also be different from zero. If D is large enough 
and negative, the system behaves in almost all the plane {T,h) as a paramagnet. 

An interesting question arises at this point: is the existence of a second tricritical 
point due to the fact that our model has a higher spin value, S = ±1, 0, than the 
Ising model, or to the fact that the S and a spin operators are alternated on the 
lattice? To address this point, we study the Blume-Capel model (BC) where a spin 
S = ±1,0 is located at each site of a square lattice. The Hamiltonian of the BC 
model is 

n = -j s,s,-DY.s^-HY.Si (20) 

<nn> i i 

where J, D and H are in energy units. We choose J = — 1 in order to compare 
with our model. Applying the Glauber-type stochastic dynamics and the mean-field 
approach described in II. b., we get that the time evolution of the BC model is given 
by, 

smh. Mm + h cos E) , , 

^ ra = -m^ — r- P ^ , ' ^ ^ (21) 

d^ cosh ^(m + /icos,^) + 2 exp(— pD) 

where m =< 5*^ > is the total magnetization. 

The numerical solution of (^Tj) shows that, when D > the BC model has an 



Ising-type behavior: it has only one dynamical tricritical point |jT2[. When D is 
very small and negative it behaves as our mixed Ising model: it exhibits a new 
first order critical line for low temperatures, as can be seen in Fig. 8. However, as 
D becomes larger but negative, the behavior of the BC model departs completely 
from the mixed model: the second order critical line disappears and, not only a new 
tricritical point does not appear, but the old one disappears, as can be seen in Fig. 
9. This study indicates that the existence of a second tricritical point is not simply 
due to the fact that our model has a higher spin operator but, to the mixing of the 
two types of spins. 

All the results shown have been obtained with f2/27r = 1.0, however we repeated 
our studies for different values of i7/27r and obtained qualitatively the same results. 
We found that the value of -D = —2.03760, at which the system starts to exhibit 
two dynamical tricritical points, is independent of the frequency, but the value of 
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D at which the model departs from the standard Ising model by the appearance 
of a second critical line depends on the frecuency as shown in Fig. 10. We found 
that, independently of the value of D, when Q — *■ the region where both solutions 
coexist (P+F) vanishes and the tricritical temperature approaches the static critical 
value. 

IV. Conclusions 

We have analyzed within a mean-field approach the kinetics of a classical mixed 
Ising ferrimagnetic model in the presence of a time- dependent oscillating external 
field. We use a Glauber-type stochastic dynamics to describe the time evolution 
of the system. We found that the behavior of the system strongly depends on the 
values of the crystal field parameter, D. For large and positive values of D the system 



behaves as the standard Ising model it has a continuous phase transition at high 
T and low h. As T decreases and h increases the transition becomes discontinuous 
and the system presents a dynamical tricritical point, see Fig. 3. However, when D 
is large enough and negative the phase diagram changes completely. First, a new 
critical line appears at low T that marks a discontinuous phase transition between 
a new coexistence region (P+F) and the (F) region. As D becomes larger (and 
negative) the ferromagnetic region shrinks, the paramagnetic region also covers the 
low h region, new critical lines appear and a second dynamical tricritical point 
emerges, see Fig. 6. The minimum value of \D\ at which the new tricritical point 
appears seems to be independent of Q (for Q 0). 

A similar study on the Blume-Capel model {S = ±1, 0) reveals that the behavior 
of our mixed ferrimagnetic model, in particular the appearance of two tricritical 
points, is intimately related to the mixing of the two types of spin operators, and 
not only due to the fact of having a higher spin operator. 

This mean-field study suggests that the mixed Ising ferrimagnetic model has an 
interesting dynamical behavior, quite different from the standard Ising model, and 
it would be worth to further explore it with more accurate techniques such as Monte 
Carlo simulations or renormalization group calculations. 
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Figure Captions 



Fig. 1. Solutions for different initial conditions with T — 0.05, h — 0.5, J7/27r = 1.0 
and D — 20. As we see, after the transient regime, some different initial conditions 
give the same solutions. 

Fig. 2. Types of solutions: antisymmetric where M oscillates around a finite value 
different from zero, and symmetric where M oscillates around zero. For T — 0.05, 
h = 0.5, Q/27r = 1.0 and D = 20. 

Fig. 3. Phase diagram in the (T,/i) plane for f2/27r = 1.0 and D = 20.0. The 
paramagnetic (P) and the ferromagnetic (F) solutions overlap in the region indicated 
by P+F. The Q symbol indicates the dynamical tricritical point. 

Fig. 4. Phase diagram in the plane {T,h) for Q/2n = 1.0 and D = —1.5. The Q 
symbol indicates the dynamical tricritical point. Notice the appearance of a new 
critical line at low T. 

Fig. 5. Behaviour of mi(0 and m2(0 for T = 0.05, h = 0.2, Q/27r = 1.0 and 
D = -1.5. 

Fig. 6. Phase diagram in the plane {T,h) for Q/27r = 1.0 and D = —2.5. In the 
inset we compare the phase diagram for two values of D: Di = —2.5 and D2 = —2.7. 
Again the Q symbols indicates the dynamical tricritical points. 

Fig. 7. Order parameter Q as function of the external magnetic field h for Q,/2t: — 
1.0 and D = -2.5; (a) for T = 0.0025 and (b) for T = 0.0305. 

Fig. 8. Phase diagram of the Blume-Capel model in the plane {T,h) for Q/27r — 1.0 
and D — —1.0. Again the Q symbol indicates the dynamical tricritical point. 

Fig. 9. Phase diagram of the Blume-Capel model in the plane {T,h) for fl/27r = 1.0 
and D — —2.5. In the inset we compare the phase diagram for two values of D: 
-2.5 and -3.0. 

Fig. 10. Dependence between the minimum value of D at which the system 
departs from the standard Ising model by the appearance of other critical line, and 
the frecuency Jl/27r. 
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